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We have developed a new on-axis digital holographic technique, heterodyne 
holography. The resolution of this technique is limited mainly by the amount 
of data recorded on two-dimensional photodetectors, i.e., the number of pixels 
and their size. We demonstrate that it is possible to increase the resolution 
linearly with the amount of recorded data by aperture synthesis as done in the 
radar technique but with an optical holographic field. 
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A synthetic aperture is used in a technique for sig- 
nal processing that combines signals acquired by a mov- 
ing detector into a unique signal field map that permits 
higher-resolution image reconstruction. This technique, 
which has been known for 30 years in the visible^ and 
the near IR^, is used extensively for synthetic-aperture 
radar— and synthetic- aperture sonar^. More recently, 
synthetic-aperture telescopes, in which several small tele- 
scopes collect optical fields that are combined to yield 
higher-resolution images, have been developed^. For all 
these syntheticaperture techniques the ultimate angular 
resolution X/D' corresponds to the equivalent aperture 
size D', i.e., to the detector displacement or to the tele- 
scopes maximum relative distance. 

Digital holography permits the complex optical field 
to be measured directly and images to be obtained by 
calculation of the field in the plane of the object, but, for 
both off-axis^ and on-axis^^^ holography, the finite num- 
ber of recorded pixels limits the resolution (which is nev- 
ertheless better in the on-axis case). To overcome this 
limitation it is natural to construct a holographic syn- 
thetic aperture by recording several different holograms 
of the same object and reconstructing the image from all 
of them. Here we use our heterodyne holography setup^ 
to illustrate this idea. Using holograms that correspond 
to various positions of the spatial filter, which selects the 
photons that fulfill the sampling condition, we performed 
an aperture synthesis in Fourier space. For a short- 
distance object, we demonstrate the synthetic-aperture 
effect by overcoming the limit on pixel-size resolution, 
which is that of one hologram^. 

Figured] is a schematic of our experiment. The signal 
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Fig. 1. Spatial filter system composed of lenses 01 and 
02, rectangular spatial filter SF, beam splitter BS and a 
CCD camera. 



field, at wavelength 1, which is scattered by a coherently 
illuminated object, passes through a spatial filter sys- 
tem, that comprises objectives 01 (/i = 50 mm) and 02 
(/2 = 25 mm), with rectangular aperture spatial filter 
(SF) in their common focal plane. 01 transforms field 
Eobj in the plane of the object into its /c-space compo- 
nents EsF in the SF plane, whereas 02 backtransforms 
EsF into EccD on the CCD. 02, SF, and the CCD, 
which can be moved, are kept aligned, and the field, 
when it reaches the CCD, propagates nearly parallel to 
the reference beam (REF). In Fig. [1] and in what fol- 
lows, the x"y" (object plane) coordinates refer to the 01 
fixed axis; xy (CCD plane) and x^y^ (SF plane) coordi- 
nates refer to the SF, CCD, and 02 moving axis, which is 
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shifted by X', Y' with respect to 01. Double, single, and 
no primes are used for the object, SF, and CCD planes, 
respectively. The resolution on the reconstructed object 
is r" = Xfi/D', where D' is the size of the SF plane 
region where Esf is measured. To avoid aliasing in the 
reconstructed images, must fulfill the sampling condi- 
tion D' < DsF = A/2/(i, which corresponds to pixel size 
d. In the best case, D' = Dsf-, resolution r" = dfi/ f2 
is equal to the size of the CCD pixels enlarged by the 
0102 expansion factor. We recorded holograms for sev- 
eral positions of the SF, CCD, and 02 moving axis and 
determined Esf over a displacement of the SF Fourier 
plane of size D' = G.Dsf (G > 1). 

For a given moving axis position {X[^Y-)^ the field 
{EsF)ij that is selected by SF is related to the measured 
CCD field {EccD)ij by 
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where z is the CCD02 distance (Fig. [T]) and Of^z is 
the lens operator that transforms field Eccd at lens 
distance z into focal-plane field Esf- The IIdsf 
dimensional gate operator accounts for the SF aperture 
and defines the size of the /c-space selected zone. These 
operators work within the moving coordinates of 02. 
Two-dimensional Dirac function S and convolution op- 
erator (g) perform the X-^Y- displacement to yield Esf 
within the fixed coordinates of 01. Operator Of^z can 
be expressed^^ as 
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X exp[i2n{xx' + yy')/Xf] E{x, y) (2) 

where Of^D Fourier transforms E{x^y)^ which is multi- 
plied by a quadratic phase function that corresponds to 
z — f propagation. To determine Esf over a wider re- 
gion it is necessary to combine the SF fields {EsF)i,j that 
correspond to various SF positions z, j: 
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If the selected SF zones are contiguous and do not over- 
lap, - = Dsf et YJ^^ -Y^ = Dsf, the wider 
region is G times wider than for a single hologram. 

In our experiment, field {EccD)i,j is measured in dis- 
crete pixels, and the CCD measured pixels are inserted 
into an X X A/" empty grid of step d. For fast Fourier 
transforms (FFTs), X is a factor of 2 larger than the 
CCD pixel number. In the SF plane, field {EsF)ij is a 
FFT calculated on a grid of step d\ Inasmuch as FFT 
steps d and d' obey Ndd' = A/2, the SF grids full size is 
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Fig. 2. Heterodyne holography setup with laser L, 
accousto- optic modulators AOMl and A0M2, beam ex- 
panders BEl and BE2, objectives 01 and 02, beam split- 
ter BS, a CCD camera, and computer PC. 



Nd^ = DsF' The IIdsj,^ S^te function is thus accounted 
for by the FFT. To obtain the sum over field {EsF)i,j 
that yields Esf [Eq. [3], we insert each N x N small 
matrix {EsF)i,j within an Esf x N' large matrix 
(TV' = GN). Accounting for the X-,Y- translation, the 
indices within the large matrix of the center of the small 
matrix z^^, j^^ are X-/d' and yj/d'. Because the object 
is at a distance z" of 01 (focal /i), we calculate Eobj by 
applying to Esf the reverse lens operator Oj^^„ whose 
expression can be deduced from Eq. [2l Because 
involves a FFT, object plane grid step d'^, which obeys 
N'd'd^' = A/i, decreases hnearly with G: d'^ = A///2G). 
In measuring CCD field Eccd that corresponds to dif- 
ferent, nonoverlapping positions of SF, one can thus cal- 
culate object plane field Eobj with a resolution and a 
pixel size that improve linearly with the amount of data 
acquired. 

Figure [2] shows our experimental setup in details. Laser 
L (633-nm HeNe laser, or 850-nm laser diode) is split 
into two beams (object and reference), which are com- 
bined by a beam splitter (BS) in the CCD camera. Two 
acousto-optic modulators (AOMl and A0M2) are used 
for shifting the reference beam by 5f = 6.25 Hz (25% 
of the CCD image frequency). A frame grabber and a 
Pentium II computer record the CCD-modulated inter- 
ference patterns and calculate complex field Ecc d in the 
CCD. A step motors (0.1-/im step) allows 02, the BS, 
the SF aperture {Dsf = 1.84, 1.9 mm in the x and y 
directions) and the CCD, keeping them in alignment. 

The first object studied is a U.S. Air Force test tar- 
get lightened by a static speckle pattern emerging from 
a ground-glass plate illuminated by the HeNe laser. In 
the CCD plane, the calculation step dx = 8.42/im and 
dy = 8.3/im is equal to the pixel size. The matrix dimen- 
sion {N X X, with N = 1024) is larger than the CCD 
dimension 768 x 576. In the SF plane, the grid step is 
d'^ = 1.84/im and dy = 1.86/im. The CCD02 distance is 
z = f2 -\- 6.87 cm, and the object 01 distance is z^^ ^ 5 
cm. To improve the resolution by a factor G = 4 in the 
x" direction we use a synthetic-aperture grid dimension 
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Fig. 3. (a), (b), the center of the image of a U.S. Air Force 
test target obtained without and with a synthetic aper- 
ture, respectively. The reconstruction was performed 
with 1024 X 1024 and 4096 x 1024 pixels, respectively; it 
corresponds to a 17.2 mm x 17 mm image. The viewed 
zone is 4.03 mm x 3.98 mm. 



of X TV, with = 4N = 4096. In the object plane, 
the grid step is then four times smaller in the x" direc- 
tion ((i^' = dx/4: = 4.21/im) and is unchanged in the y" 
direction {dy = 16.6/im). 

We calculated for a single hologram the U.S. Air Forcer 
intensity image (N x N = 1024 x 1024 pixels, or 17.2 mm 



X 17 mm). The central part of the image is shown in Fig. 
[3l^a). The image resolution and the speckle size, which 
depend on /c-space extension = Dsf , should both 
equal d" = {fi/ f2)d. In our case, a low-pass filter, in- 
ternal to the camera, lowers the bandwidth of the analog 
video stream of data, which degrades the resolution and 
enlarges the speckle size by a factor of ~ 2 in the x'' 
direction. 

When the synthetic-aperture algorithm is applied, be- 
cause the BS moves with SF the reference optical path 
length changes for each SF position. Each measured 
EccD field map is thus shifted by an unknown phase. 
To determine the phase correction we acquire holograms 
such that the zones covered by the aperture overlap 
for two consecutive SF positions {X[j^^ ~ ^i) < DSF 
{i = 1...21 and X-_^-^ — X- = 250/im in our experi- 
ment). We assume thus that fields {EsF)i-\-i,j are equal 
in the overlapping region, and we determine the phase 
by maximizing, in the overlapping region, the {EsF)i,j to 
{EsF)i-\-i,j field correlation, which we found to be 90% 
at maximum. We calibrated the X-_^-^ — X[ transition 
(135 pixels) and the CCD02 distance [z = 6.87 mm in 
the quadratic phase factor of Eq. [2] by this optimization 
process. 

We calculated the synthetic-aperture reconstructed 
image (4096 x 1024 pixels, 17.2 mm x 17 mm) when 
the SF, 02, the BS, and the CCD are X' translated 
over D' = 5 mm in 21 250- jam steps. The image cen- 
ter is shown in Fig. ^h). The motor calibration yields 
a relative translation of 134.8 pixels, in good agreement 
with the value obtained by correlation. As expected, the 
synthetic-aperture image exhibits better resolution and 
smaller speckle size. If the BS phase-shift correction is 
not made, the speckle size, which depends on the field 
extension in /c-space ~ GDsf^ is the same, but the res- 
olution, which depends on the /c-space field coherence 
length (~ GDsF with correction and ~ Dsf without), 
is lower. The image contrast is thus lower. 

To perform a quantitative analysis of our synthetic- 
aperture technique, we studied a narrow hole of diameter 
= 30 mm located at z = 5.5 cm in front of the SF. 
The hole is backilluminated by an 850-nm laser diode. 
The SF, 02, the BS, and the CCD are X' translated 
over D' — b mm in 250-/im steps that correspond to 
101.6 pixels. 01 is removed, and two Fourier transforms^^ 
propagate the field from the SF plane to the object plane 
such that the grid size ( d'^ = d'^ = 2.46/im, dy = dy = 
2.5/im is the same in both planes. Figure H] presents, in 
the x^^ and y^^ directions, the reconstructed intensities 
(points) compared with the theoretical intensities (solid 
curve) that account for pixel averaging and diffraction. 
As expected, the edges are sharper in the x^^ direction. 
For both x^^ and y^\ the agreement with the experimental 
points is excellent. This experiment, performed without 
01, can be interpreted as yielding the synthetic aperture 
in real space for field Esf that we obtained by Fourier 
transforming ECCD. The synthetic aperture here allows 
the images angular resolution to be increased to reach 
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Fig. 4. Cut in intensity for a 30-mm hole in the x and y 
directions. Points, experimental synthetic-aperture re- 
construction; and curves, expected theoretical shapes. 
Dark filled squares and the solid curve correspond to the 
x" cut; lighter filled circles and the dashed curve corre- 
spond to the y" cut (synthetic aperture is formed in the 
x" direction). 



grid step limit dO for the resolution on the image. One 
could also use the real-space synthetic aperture in the 
CCD planell without transforming Eccd into Esf^ but 
the ultimate resolution limit d is much less d' <C d. 

The experiments presented here are examples of ways 
in which synthetic apertures are obtained in both real 
and Fourier space. In Fourier space, the synthesis allows 
either the field of view in the far field^ (with respect 
to Nd^ /\) or the resolution for the near field (U.S. Air 
Force experiment) to be improved. In real space, the 
synthesis improves either the field of view (near field) or 
the resolution (far field). In the experiment with holes, 
the pixel size of field Esf-, d' ^ is small, and the hole is 
pushed into the far field z" > Nd''^\. Note that for both 
the U.S. Air Force experiment of Binet et al.^^ and the 
hole experiment, the synthesis yields better resolution, 
not the enlargement of the field of view that could be 
obtained by a simple scanning method. 
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